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ABSTRACT 


This report investigates subsonic diffuser performance, with 
emphasis on conical and annular geometries. A correlation is 
presented which aids in the prediction of performance. Two an- 
nular diffusers were designed and tested to help substantiate the 


correlation. 
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I. INTRODUCTION 


The diffuser is a device used in fluid mechanical systems to 
convert kinetic energy into static pressure by decelerating the 
flow. There are several basic straight-walled diffuser shapes 
such as conical, two dimensional and annular as shown in Figure 1. 
Much attention has been given to two dimensional diffusers and 
several successful correlations have been developed, but annular 
and conical diffusers have not been as exhaustively studied and 


good new correlations would be very useful to the designer. 


A. A GENERAL DESCRIPTION OF THE FLOW 

A diffuser should convert kinetic energy of the flow into 
pressure with a minimum of lesses due to viscous effects. Also it 
is desirable to make optimum use of the geometric area ratio, 
because the area ratio prescribes the amount of diffusion that is 
possible under ideal conditions. Bernoulli's equation relates 
the total fluid energy to the static pressure and kinetic energy; 
and if the flow experiences no losses due to friction or other 
non isentropic processes the total energy of the fluid remains 
constant. Any decrease in velocity will cause the static pressure 
to rise, so that in an ideal diffuser all kinetic energy losses 
show up as pressure gains. Unfortunately, most fluids experience 
losses, and the goal of any flow analysis is to keep track of the 
fluid energy whether it is dissipated or converted from flow 


energy into static pressure. Dissipation of kinetic energy takes 





place inside the boundary layers where large velocity gradients 
exist so that the viscous effects become important. 

The boundary layer is a thin region adjacent to the wall in 
which the velocity increases from zero to the free stream value. 
Since the flJuid in the Bente layer is moving slower than that 
in the free stream, it does not have as much kinetic energy as the 
bulk of the flow. The transverse velocity gradients that are 
present give rise to shear stresses which are related to the veloc- 
ity gradients through the coefficient of viscosity for laminar 
flow. In a turbulent boundary layer the shear stresses are not 
related to the gradients in Such a simple fashion. However, it is 
possible to model turbulent flow as an average Tee earelus time 
dependent fluctuations whose average value is zero. When the 
boundary layer equations care modificd for turbulent flowe, addi- 
tional shear stress terms are present known as Reynolds stresses. 
These Reynoids stresses are in fact the predominant stresses for u 
turbulent boundary layer. 

Since the flow in diffusers is constantly decelerating, due to 
increasing area, the flow faces an increasing pressure as it 
traverses the channel. To overcome the increasing pressure the 
flow must transform kinetic energy into pressure. However, as 
mentioned the flow in the boundary layer is not as energetic as 
the free stream, but it feels the same longitudinal pressure gradi- 
ent. When the fluid near the wall does not have the energy to 
overcome the pressure rise, the flow higher up in the boundary 
layer must pull it along by shearing action. In this sense the 


turbulent boundary layer is more capable of meeting an adverse 


LO 





pressure gradient than a laminar one because of the additional 
Reynolds stresses. For this reason it is advantageous for dif- 
fusing flows to be turbulent in order to withstand the decelera- 
tions imposed in a diffuser. A point is reached, however; where 
the outer layers can no longer pull the inner layers, and the flow 
is forced away from the wall and separation results. Large scale 
separation with recirculating flow causes intolerable losses in a 
diffuser and is to be avoided. Any time the boundary layer does 
not have sufficient energy to overcome the presSure gradient it 
will separate, so it is necessary to find out something about the 
way the pressure gradient behaves in a diverging channel. 


Bernoulli's equation is written as follows 


R= ree +t (1) 


Since the boundary layer feeis the same pressure Gistribution as 
the free stream, it is only necessary to calculate the free strean 
pressure distribution which shall be calculated here assuming 
inviscid one dimensional flow in the core. Taking the derivitive 


of (1) gives 

dp/dx = - vdV/dx ™ 

with continuity ——— density 

mw“ = fAV = 4, AV, | (3) 
or 


dA/Az -d¥/V om vay 


) 


(-v*/A\ oA ") 
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replacing (4) in (2) 
dp/dx = (-7V"/A) dA/dx (5) 


Fech COmelnuicy 


dp/dx zs (7 A," V,*/ A?) dA /dx (6) 


ACHE pyD/dx = ACH) fdx = (2AS/a3) dA/dx 


(7) shows that the pressure gradient is the strongest at the inlet 
of a diffuser since the area gradient, dA /dx , is very nearly 
constant for a given straight walled diffuser. This is compatible 
with the state of the boundary layer which is “youngs" at the en- 
trance of the diffuser and is able to overcome large pressure gra- 
dients. The large initial gradient suggests that a boundary layer 
Should be as thin as possible at the entrance to a difiusetr tor 
best performance; i.e., avoidance of early separation. This fact 
has been substantiated by many researchers by increasing the inlet 
boundary layer thickness and noting a decrease in performance, 
Kline [Ref. if Returning to the expression for the preSSure gra- 
dient, (7), it can be seen that the gradient may be eh oni btariys 
But a large area gradient also insures that the pressure gradient 
term will decrease rapidly, because the denominator of (7) contains 
the local area raised to the third power. Although a general 
description of the flow conditions in a diffuser has been given, 
the task remains to determine losses due to the factors considered 


above. 


le 
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B. REVIEW OF PREVIOUS WORK 

Many attempts over the years have been made to explain and 
predict the actions of the flow in the diffuser. As the previous 
section suggests, it is not a simple process. The major stumbling 
block is that the phenomenon of turbulence is not well understood, 
and to predict separation the growth of the boundary layer must be 
calculated. In order to bypass the turbulent flow calculations, 
more empirical methods and simplifications have been used to pre- 
dict performance. One of the simpler techniques is that of the 
equivalent cone angle method as described in Gleason [Reg. ul. 
Basically, all flow cross sections are related to an equivalent 
conical flow eh a corresponding wall divergence angle. If the 
wall divergence angle, which is a direct function of pressure gra- 
dient, is below a certain value the diffuser is judged sound. 
This method is widely used for design but has not always proved 
adequate. The next level of effort has been directed towards gen- 
erating performance plots. | 

After testing a sufficient number of diffusers of a given type, 
constant pressure recovery lines were Se eted on a graph of perti- 
nent geometric parameters such as wall angle versus length. 
Kline Vref. 5] has done extenSive work in eee for two dimen- 
sional diffusers and van and Klomp [Ref. 6 | have done the same 
thing for annular diffusers. Through flow visualization techniques 
Kline [Res 5 | has been able to generate a malo showing the expected 
elle regimes for different two dimensional geometries, Figure 6. 
The line of first stall, a-a, is presented as the optimum perform- 


ance line, because at or just after this point the boundary Laye? 


i 





has overcome as much of the pressure gradient as possible and is 
nearing large scale separation or transitory stall. No similar 
flow regime graphs have been developed for the annular or conical 
geometries | 

Several years ago Kline [ Ref. 1] developed a more analytic 
approach to the problem. He has attempted to solve the turbulent 
equations in predicting performance. The boundary layer and con- 
tanuity equations have been combined into a mementum integral 
equation and used with several shear correlations to work up a set 
of five linear first order differential equations with non linear 
coefficients. However, the solutions of these equations alone do 
not predict the performance, a stall criterion is needed. Kline 
has combined several parameters including boundary layer thickness 
and rate of area growth into a parameter. When this parameter 
decreases to a certain vaiue Siaii is said to occur. The param- 
eter has no physical basis, but it does correlate line a-a of the 
flow regime chart. In Reference 1 many reports have been analyzed 
and the criterion in fairly successful for 2D-diffusers. The 
variation in performance with inlet boundary layer thickness -is 
also demonstrated. Though Kline appears to be fairly successful, 
the method is rather approximative and requires the solution of a 
large set of equations. However, Vavra (Ref. 8 | has developed a 
performance parameter which is easy to calculate from geometry and 


can be derived from basic principles. 
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Ii. DERIVATION OF THE SHAPE PARAMETER OMEGA 


Vavra (Ref. 8} presented the basic derivation, and Gapp (Ref. 3] 


did a thorough development which is presented below. 


The momentum theorum 


JL din V- f dm, Va -f@pds- fiends (8) 
Gy) 2 & 6) ed 5 5 


was applied to the differential flow element of Figure 2 to yield 


(fav) Vedv)- Gav) V = 


= & oe A+ da)on, pdA -mpdd- Foc ( 
2, pA a, (prap) ¢ oe ® 
After expanding and rearranging terms where Ay A fy Gnd “My, are =n 


oun» 


outward pointing normals of each respective side and & is the 


tangent vector, the theorum gives 


fAvdV = ~Adp- 2cdl (10) 


rearranging 


dp = -pVdV- ~c dl /Aa and 


or 


G (Vv*/2) e dp /+ = (- Z/¢\ CCA) dl. (1.2) 


Equation (12) may be derived for an annular diffuser as well. 


From Figure 2 (a) for an annulus 


dL,= dh cosy ) (13) 
A, = AfSCOS ¥ (14) 





Prom Figure 2 (b) 
ae eir CF, +2). = Heat 4 (15) 


Also Cl! =4r Pq tO¥ small divergences. Replacing (13), (14), and 
b>) an (12) 


d (vipa) + dpyy = (- @/4) (C /A,CoS &) dk, cost = 
wre! ah, (PA, 


So finally considering only sections normal to the axis and 


omitting the prime 


d (V'/2) +dpsf = C-T/S)CCfA) dh . 


as before. From the first law of thermodynamics 


oS (16) 


~. 


° ® 
oo a on af ff 


. f _. 
6Q= Gut par -yap= dh- 7 op= 


- ah 
P 


. ee 


or 
dp / 4 = dk - Td (17) 
Replacing (17) in (12) gives 


dh+ d(V*/z)-Tds = (-2/7) (esa) 4b (18) 


combining the enthalpy and velocity derivatives into the total 


enthalpy form 
dt = Tds- (€&/f) (C/A) db (19) 


Since the idealized one dimensional process is adiabatic, with no 
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energy input through the walls, the total enthalpy is constant and 


its derivative is zero. Therefore (19) reduces to 
Tds = (t/p) (C/A) dL. (20) 


The shear losses in a turbulent boundary layer are not easily 
determined, but they do give an indication of overall losses. An 
average shear stress coefficient was defined for the entire chan- 


nel to reflect these losses. 


Co = Cav Si 7 V~ (21) 


Srsusing continuity 


Say = (Cg/2) m*/f a” (22) 


Replacing (22) in (21) yields 
a 2 née 3 
Tds = (Cg /2) Cn LFA )¢ dk. (23) 


Now a Similar differential entropy change will be developed from 
purely thermodynamic considerations for an adiabatic polytropic 

compression. Referring to Figure 3, the polytruopic efficiency is 
defined to be “Fp = dT}<or uSing the isentropic relations for a 


ay 


differential compression ¥ [ 
mT. = ve (24 
(T+ dTj.)/T ro et dp/p) — 
Expanding in a binomial series and neglecting higher order terms 


“feet 
(1+ spp) =F + dp, hugher once renme (es 


Replacing (25) in (24) and rearranging 


dtssdt = (T/dT) psp) (26) 
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using Yo = dTis/dT 
dT/T = Cer) yp) doy p (27) 


Introducing the perfect gas law and p= pe into the first law 


of thermodynamics yields - 


Tds = RO¥/e- dT (RT/pP) dp 
@s= R { ¢ ¥/y-) ad Tyr - psp] (28) 


The final differential entropy form results by using (27) 
ma (eo) 
=, t 
dye = (Yep -1) dpyp 0 
In order to combine the geometric and thermodynamic developments 


the differential entropy changes, (29) and (23), can be equated 


(Cp /2T) va® Cdl =  R(My2-I) dpyp OD) 
T*A roi 


. ona a 6 
Introducing the dimensionless referred mass flow rate > io ae 


into (30) yields Pi, 


(Co/ 27) ( ya A= ee i/ RT, 4* A: ) C wk. a 
(t/ 7) ( (> -1) ap / p (31) 

again using the perfect gas law f a ww RT 
(Pip) Ce 27) Cre APA) Cd Lb = 
(T/f) ( Yap “Op (32) 
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rearranging 

(C72) (bay Ar /A3) CAL = (PA) (Mp-d(P/p) (33) 

For a polytropic process Fyn = constant so (33) — 
(G/2)m> (A>/A*) Cdb = C4 p -1)¢/p,)" “d Pip) 4 


The equation is now in a suitable form to be integrated. Vavra 
has called the left side dX, and integrated it from the channel 
entrance to exit along s streamline so that the constant total 


eye PY go DETOR is valid. 


f dx, = [C%/2) m2 (A,*/at) Cdl @s) 


As previously .discussed C. is a constant representing, shear losses 
and may be brought outside the integral. Also referred flow is a 


constant for a given flow rate. a 
m h 


we ge fn ; ae ff ‘Se fun 2 , 
J ah, CCg /2.) m, JX 


t 4° 


«e* , 


on iF t SUN 
\eo i» (YU) 


Nee © 


The integral of (36) contains geometric characteristics of the 


diffuser and bas been defined as the shape parameter omega (aJis) . 


Qe ¥ ‘Athy Cor AS db. (37) 


The product £}, Ce must be related to diffuser performance for 
a). to be useful. In order to do this the right side of (34) 


must also be integrated. 
AX, = ( Vepnt) C Pp,) ™ dC P7P,) 


Integrating as before from entrance to exit. 


F/P, fee 
x, f "OY -1) C ) d CF7?,) (38) 


ii 


Pall 





By assuming constant polytropic efficiency and n constant 


W/o 
x, = (4% p -1 ) Ch/ ni) (cPupy” ~1) (39) 


For the differential change Po = Pit oP with of é¢ Pp, and 
again using a binomial expansion 


net 
[ (P+ Ap) jp, ] = = b+ Cn+rif;) AP /P. + higher order terms 
replacing in (39) 
kK, = (Yay —!) Psp, (40) 
equating (36) and (40) 


b -2 
COgja) vane [(Ava) Cla) dh = 


CF /2) vw n* DL = C in » —1) QP/P, (41) 


Substituting Tor mr,according tu tie definition 


IM, » “fi AV, VRT, /P,A, 


gives 


(C, 2 CF, AP V> RT) 2 /0P,* A) = Cp a AP/P, 


Na Ce = CLAP/ FM") C Yap -12 


Since constant total enthalpy was assumed there is no change in 


(42) 


total temperature so rewriting (28) 
dg = -R (dR /) (28a) 
equating (29) and (28a) 

-dR/P, = CV 1) epsp 


(29a) 
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(29a) is in a suitable form to be integrated and is a straight- 


forward a“ P./P. 
(- vn Pe), em ( (1 —}) i | Pe 
In ¢ Pa, / Pa, ) = ¢ (ao-l) li ei. {ep (43) 


again using 


Pe J. = A, /CP,4 OP) ands Pap, = It APZP, 


with the usual expansion for logs; i.e. In(l1+x) = x, (43) becomes 
i ee = = 

App, ( ’ep-/) APs / Px, (44) 

or 


CYap-!) Ap = ~ Op, CPi /m,) 7 AR (Oe Em) 9 


All flow considered was incompressible with M<<L SO 


(Ya, “')Ap= —aP, (46) 
replacing (#6) in (42) 

D.Cg = -2 Op CFV) = CP a ~ PL YCE AY 

Pe Pay CEs 2) = CP ea Bette WELW) 


(48) 


Me 
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However Cp = coefficient of pressure 
Cp = CPP) / Ce (49) 


E>: = pressure recovery if the process followed the P, to 


G 
P line of Figure 3, i.e. ideal. 


iSen 
Cp. = C? =-BY/ GS Y%*) (50) 


om lO 


For isentropic flow AP, =O so (50) becomes 


Cp. = CF/2M," ~ Sia vo egy) = Io (Ai/A,)* (S02) 


Replacing (50a) in (37) gives 


CPe, Pe V/C ELV) = 


ab ye | a z 1/6 *) 
using (49) and (48) in (51) yieids 
= Coy. ~ (52) 
ILCy Pe Cp 
(52) is the basic result of the theory. The amount a flow 
departs from isentropic conditions has been shown to be a product 


of the geometric parameter omega and an average shear stress 


EoernreGient . 
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FIGURE 3 


I-s diagram foran adiabatic compress ion 
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III. ORIGINAL SUGGESTION FOR USE OF OMEGA 


Once equation (52) was derived, Vavra (Ref. g ] investigated 
some two diiensional data by plotting the losses as a function of 
the shape parameter. His plot is reproduced in Figure 4. It can 
be seen that the losses tend towards a minimum value in the neigh- 
borhood of an omega of ten. For this reason Vavra suggested that 
an omega of ten might be a good diffuser design value. 

In an earlier investigation of the shape parameter, 

Gapp (Ref. 3] reduced the: annular data of Sovran and Klomp [Ref. of 
so that a plot similar to Figure 4 could be generated as shown in 
Figure 5. Although the plot shows the same trends as Vavra's 

' graph, there does not seem to be a value of omega to insure mini- 
mum losses. Additionai two dimensional data of Kiine (Ref, 5} was 
also plotted and has very similar characteristics, but does not 
offer any firm evidence that one value of omega is to be strongly 
preferred a any other. A somewhat different approach was to 
plot lines of constant omega on the two dimensional flow regime 
chart Fox and Kline (Ref. of, Figure 6. It was seen, however, that 
a given value of omega did not necessarily indicate the state of 
the flow. Although the above evidence indicates that Vavra's 
original suggestion of how best to represent losses as a function 
of omega was not the best implementation of the shape parameter, it 
does not adversely reflect on the validity of omega as a good 
parameter. Omega has a definite physical significance which is 
Fundamentally different from previously suggested parameters 


embodying all pertinent geometric characteristics. A der. sation 
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of the physical meaning of omega has already bcen given and an 
additional aspect of omega will be given in a later section. 

The effort to achieve small losses as reflected in a minimum 
value of tts Cy may not always be the chief goal of a particular 
diffuser. The ideal pressure recovery coefficient is a function of 
area ratio and as such is a measure a the maximum diffusion ob- 
tainable for a particulr geometry. If the designer wishes a dif- 
fuser to convert a large portion of the flow evergy into pressure, 
a larger arca ratio is called for. However if only a modest rise 
is called for with a fairly uniform exit profile the reverse is 
true. The application calling for a large pressure rise may be 
able to tolerate higher losses when maximizing recovery while the 
need for a uniform velocity profile will tend to require minimum 
losses In pronosings a theory or correlation for performance the 
two above conditions must be considered so that an optimum design 
criterion may be formulated. Before investigating diffuser 
characteristics further a somewhat simpler derivation of omega will 
be given to relate the shape parameter to diffuscr performance in 


a different way than previously derived. 
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IV. AN ALTERNATE DERIVATION OF OMEGA 


The non dimensional group used to define omega was a logical 
consequence of the combination of the momentum theorem for an 
arbitrary channel and the expression for a differential entropy 
rise in an adiabatic polytropic compression. The derivation was 
necessarily lengthy and the expression for omega may be obtained 


from a more baSic approach. 


For the same one dimensional channel flow of Figure 2 


oo fe tf Hp (53) 
n= SAV = SAY, su) 


Taking the derivative of (53) with respect to iengtn and assuming 
constant total pressure 


dP, = dp +. pVdVF0 (85) 


are@eiecom continuity for constant density 
VdA = -~AdV 


or 


vdV = (-V*/A\ d A (56) 
replacing (56) in (55) gives 


% =, (37) 
fV" dasa =dp 
and using the definition of mass flow from (54) 


G ke 


a 





rearranging 
dp/ cd pyr = 2 CA,/A\* CfA) CoA) dL 6%) 


and integrating both sides 


Cp. = 2 fA, 7A) (YA) CAAL) dL (60) 


L 


The integral in (60) looks very much like the definition of omega, 
but the area rate of change must be related to the wetted perimeter 


for each diffuser type. 


A. CONICAL 
a 
Achy = 2m (R, + LTAW®) (61) 
dA/d. = 2r TAWve CR,+ LTA) (62) 
Coby = 277 CR et TANS) (63) 
PiGinine (63) and (62) yields 
C(.) = (dA/dt) /TAYE (614) 


and replacing (64) in (60) 
‘ a 
CP, = pat , (A (JAAD COL) TAWVEO dh =2THVO 2 (65) 


Therefore for the conical case the ideal pressure rise is a 
function of the wall angle and omega. The factor relating 


CPi to omega has been defined as PARHC. 





B. TWO DIMENSIONAL 


dA/dL = 2B, TAWE (66) 
C(t) = 28, + 2w,+ 4 LTAV6 = 2B, + 2welly) 


2B, = ccéLy) -2rawee) (67) 


combining (66) and (67) 


dA/dL = CCL) TAVOE CJ] ~2Wwlt)fecr)) (68) 


Cp; — {rrawe Ci-2wees/eer)) (& i *( ) sh com 


Equation (69) is also similar to the expression for omega, but an 


investigation has shown that Sp. is not related to omega by a 


Simple constant analogous to that of the conical diffuser. 


C. ANNULAR 

An annulus is a complicated shape, but it shares certain © 
Fundamental. geometric characteristics of the cone such as axial 
symmetry. In fact, an annular diffuser is nothing more than a 
conical diffuser with an Belsyvinne bile cen tLenebody .alOrsanmannuian 


diffuser: 


dA/di = rh (CTAV> @, ~TAWV* 6c) 
+ 297 (RT) TAVO, - RHI TAV OE) (70) 


CétL) = 27 ( Bett) t Rec] ae 





dAasdi = [C Taw ©, RL) -TAW 6, Ry ()) Cau) (72) 
CRelty + Rye CL) 


By examining the factor relating area gradient to the local cir- 
cumference in (72), it can be seen that an expression similar to 


(65) can be derived for annular diffusers. The factor is defined 


as PARHUA. 


PARUA= [TANG Rolt) -TAVO; Rr(U]/C Rete) +2Q)] 


By referring to Figure 7 and making use of the theorum relating 
Similar triangles from geometry it can be seen that PARYA may also 


be expressed as a function of the inlet geometry. 


~ A Aa —_ SA a i od ” t <s 3B ° bid 
PARYA = TAWE, (CATI/Lin} Bo TeV OS CBHI LL )L 


Cy se ewe 


CR, / ui L. + (R Thf Lis) t. 





wwe To et 


tri 


= (TAVO, -RRAT TAVEeL) (/CIt RRAT) 


RRAT = inlet radius ratio = QQ. He / Lo 


It can be seen that PARYA is constant for a given annular diffuser 
family. A family is defined as a set of diffusers with constant inner 
and outer wall angles and the same radius ratio. For example all 
diffusers with wall angles of 15 and 20 degrees with a radius 
ratio of one half would comprise one family while all diffusers 
with wall angles of 10 and 25 degrees with a radius ratio of ./ 
would comprise another family. The previous derivation may be 
used in (52) 

SL. Cy - Ca, ~op (52) 
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or 


WC, = Cpe C1 Sp/cpi) 


Co = (Cp. srs) Ci- CP f/ep,) = 
PARY Cl - CP/C_p;) (52a) 


Although the previous derivation of omega does not allow direct 
Calculation of performance, it does give a better understanding of 
the shape parameter. Once the above derivations were formulated 
a parameter study was undertaken to observe the eeenonce of omega 


to certain geometric factors for each diffuser shape. 
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(INLET SECTION 





FIGURE a 
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V. PARAMETER STUDY 


The parameter study was done with the aid of several computer 
programs which are listed and explained in Appendix A. To doa 
complete study would have been a major undertaking with all the 
associated variables of conical, two dimensional, and annular dif- 
fusers. However, a representative study was done by varying each 
parameter a reasonable number of times to establish trends. The 
conical diffuser was the simplest shape and was studied first. 

The conical variables were taken as inlet area, length, and 
wall angle; but to have a manageable amount of data the inlet area 
was set equal to unity and the other parameters were varied. 
Figure 8 shows that omega is a strong function of wall angle and 
past a certain non dimensional iengtih, the wail angie dominatcs. 
For a high divergence angle the characteristic length is short be- 
fore a constant omega is achieved while the opposite is true for 
the lower divergence angles. The wall angle is inversely propor- 
tional to the final value of omega as shown in Figure 8. An 
interesting fact about Figure 8 is that an omega of ten results 
‘for a wall aie Cre avOl ithe cm@euneces = Cenredte da tr Use tS) tier 
small divergence angles have low losses, but they do not neces- 
sarily have the highest pressure recovery for lengths common in 
practice. A better look at omega in relation to flow conditions 
was found by plotting constant omega lines on the two dimensional 
flow regime chart of Fox and Kline (Ref. 5]. 

Figure 6 was briefly mentioned earlier, but there are several 


aspects of the plot that were not discussed. [Tor the two 
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dimensional geometry, omega was found to be dependent on angle, 
length, and inlet geometry. As Figure 6 demonstrates, a particu- 
lar value of omega in the two dimensional regime can indicate a 
full range of flow conditions from unstalled to large transitory 
Seat. In addition to Pe indicating conditions in the two 
dimensional diffuser, the constant omega line takes on different 
values depending upon the inlet geometry. 

Finally the characteristics of annular diffusers were studied. 
The annular diffuser is the shape of most interest for turbo- 
machinery applications and unfortunately is the least understood. 
There are numerous parameters to vary: wall length, radius dif- 
ference, mean radius, inner wall angle, and outer wall angle. 
Due to the large numbcr of variables no plots analogous to 
Picure 6 and Figure 7 were apparent to show omega trends. However, 
one graph which combined inlet geometry and wall angles was plot- 
ted versus non dimensional length. Figure 9 has the advantage of 
incorporating as many geometric parameters as possible into a 
meaningful group by plotting area rate of change versus length at 
a constant value of omega. The area rate of change is Sane come 
stant for a given set of wall angles and inlet geometry as can be 
seen from (70) when € TA nye = -TAW*6}can be neglected, 
and it is also related to the pressure gradient. The plot shows 
the same general trends as the other two geometries by indicating 
the omega dependence on area rate Bp orenenes and length. The 
figures all show that high omegas can only be achieved with low 
area rate of change, but several geometric conditions may result 


~~ 


in low values. A diffuser may be very short or have a hi‘h are 
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gradient in the case of a low omega value. Although the parameter 
study indicated geometric trends, it did not offer any insights 
into performance. To gain an understanding of diffuser perforn- 


ance the published data was analyzed more closely. 
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VI. DATA CORRELATIONS 


Previous attempts to predict or explain diffuser performance 
relied on a number of parameters and graphs that were rather 
unweildy to use, or an oversimplified approach such as the equiv- 
alent cone angle method was employed. Gapp (Ref. 3] looked at the 
data of Sovran and Klomp fRef. 6] in the manner suggested by 
Vavra [Ref. 8] and tested several diffusers, but he was unable to 
develop a method of prediction based on his work. It has been 
shown in previous sections that the parameter omega has physical 
Significance in relation to diffusing flows, and it was reasonable 
to assume that omega is in some way reJ.ated to performance. If 


some means of determining the shear term, Cr, could be found; 


<! 
1) 
rs 
C) 
[- 
9) 
nr 
cr 
() 
3 
2 
al 


equation (52) would vive performance directly. Sev 
were made to arrive at a technique for determining Ces but some 
other methods of predicting performance from curves bear mention- 
Pio 1 Erst. 
If ~¥.Cp is small the recovery is very close to ideal, but 

aftwCp must be minimized for a given area ratio es nerd maximum 
pressure recovery. This fact suggested a possible utilization of 
the original plot of Vavra. By drawing lines of constant area 
ratio on an-.ft. versus lh Cr graph, an optimum omega range can 
be determined at the given area ratio. For the given area ratio, 
the omega range would guarantee minimum losses and maximum pres- 
sure recovery close to the desired one. A plot such as Figure 10, 


which is a line of area ratio two from the data of Sovran and 
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Klomp [ Ref. 6), would result. It shows an omega range of five to 
eight for minimum losses. Unfortunately Sovran and Klomp did not 
have a wide range of area ratios, and to do the required testing 
to generate other data would be very costly and time consuming. 
Another plot was discowenad, but it did not have a wide 
range of applicability. However, it is worth mentioning since it 
correlated a large amount of the annular data of Sovran and Klomp. 
The non dimensional area gradient, ( dA/dL) JAR, is related to 
the pressure gradient, as shown in (57). After several attempts 
to relate the ree sradient to losses a graph such as Figure ll 
was plotted. It is a plot of non dimensional area gradient versus 
losses for a given non dimensional length, and a Sees line 
resulted for a range of lengths. The resultant straight lines 
Sm@meeested clcts such as Figures 12 and 13 in an attempt to cor- 
relate all the annular data. Figure 12(a) is a list of symbols 
awa in Figures 12 - 15. However, considerable scatter is evident 
in addition to some definite trends. Two distinct lines are vis- 
ible and some data considerably off both curves. The two distinct 
lines correspond to families with inlet radius ratios of .55 and 
.70 respectively. The data below the lines are from families with 
inner and outer wall angles that differ by no fore than two degrees. 
These plots are of limited value at best and are rendered even 
more questionable because the diffusers with the best recovery, 
the families with similar inner and outer wall angles, do not even 


fall on the curves. 
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The independent derivation of omega introduced a geometric 
parameter that related the ideal recovery to omega (see 
equations (52) eee (72). The parameter, PARY, had a different 
form for each class of diffuser, but in all cases it corresponded 
to the quotient of the area gradient and circumference for a 
given axial station. If the shear stress coefficient could be 
determined, the theory of Vavra would give performance directly 
from equation (52). To see if the parameter, which has been desig- 
nated PARYA and PARHYC, had any relation to the shear stress coef- 
Ficient, Figure 14 was plotted with the conical data of Cockrell 
and Markland fRef. 2} and annualar data from Sovran and Klomp. 
The conical data consisted of twenty-four diffusers which were 
tested at two boundary layer thicknesses and approximately eighty 
annular diffusers. No information was given concerning the exact 
dimensions of the inlet boundary layer for the annular tests, but 
femetbe annular diffusers were tested on oe same apparatus with a 
thin inlet boundary layer. The set of conical diffusers which is 
displaced upwards on the curve had a thicker inlet boundary layer 
and exhibited the expected higher losses. Some scatter is Syitene 
in the figure for the longer annular diffusers, but a definite 
curve is discernable. 

Once Figure 14 was established, the next logical step was 
Figure 15. To enable the direct determination of recovery from 
a known omega, PAR4 was plotted varoue the quotient of pressure 
recovery and omega. Again a good curve with small scatter resulted. 
The conical data with the thick inlet boundary layers were dis- 


placed downwards as expected. Tigure 15 differs from the previous 
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annular correlation in that the data fall close to the line, and 
the wide scatter attributed to families with similar inner and 
outer wall angles is not present. Also the curve does not seem to 
depend on inlet radius ratio. With Figure 15 and a knowledge of 
geometry for a conical or annular diffuser with a thin inlet 
boundary layer thickness, the eee recovery can be determined. 

Besides predicting the recovery, the correlation presented in 
Figure 14 can be utilized in the development of optimization 
criteria. The designer is often faced with the problem of maxi- 
mizing the recovery for a given length. The annular case will be 
discussed, but the development for a conical diffuser is similar. 
It was shown that 


LLC = Co- 


G 


Combining the equations gives 


Cp = Cp. Ci- Ce /par4) (73) 


Q 


However, Figure 14 shows that C- is a function of PARY. The curve 
is of the form Cr -? Kx OPA p. 4) : , Where cf ead b&b are 
constants that can be determined from polynomial curve fitting 
techniques used in numerical analysis. Also the ideal recovery 
coefficient can be shown to be a function of inlet conditions and 


non dimensional length. 


Cp, a | — (4.,/A,)™ 


~~ 
AJA, = CeTi+RHi)/Lamvr+2as rae C18 | 
CTé# ef “0, = Tid “9, ) t "dt bes ( oe PT Oe aM 


(74) 
RHI TAV OE; > | 
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By replacing (74) in the expression for the geal recovery coef- 
ficient and combining this result with an expression for Cr as 
suggested, the pressure coefficient will be expressed as a function 
of inlet geometry and non dimensional length. When an expression 
has been developed, the theory of maxima and minima of a function 
may be used. By taking the derivative of ©s with respect to non 
dimensional length and setting the resultant expression equal to 
zero, the recovery could be maximized. To obtain a given recovery 
in a minimum length, the correlation could be employed in a similar 


fashion. 
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VII. EXPERIMENTAL WORK 


A. EXPERIMENTAL EQUIPMENT 7 

Two annular diffusers- were designed and tested in order to 
provide further data in the evaluation of any trends or correla- 
tions discovered from analysis of the literature. The diffusers 
were designed with the aid of program DESIGN which is listed in 
the appendix. With a given inner wall angle, fifteen degrees, the 
program varied the outer wall angle and length to maintain a given . 
area ratio and calculated omega for each case. Both diffusers had 
an area ratio of three and an inner waJll angle of fifteen degrees. 
Model one had an axial length of 7.73 inches with an outer wall 
angle of nineteen and one half degrees, and model two had an outer 
wall angle of eighteen degrees wilh an axiai length of 9.36 inches. 
Both diffusers had an inlet hub radius, Rhl, of 1.875 inches and 
an outer inlet radius, Rtl, of 3.123 inches. Model one had an 
omega of four and model two had an omeza of five. Both were con- 
structed from phenolic resin. The above information is summarized 
in Figure 16. Figure 17 is a photograph of the a ody and 
model two. In order to supply a uniform flow to the diffusers, a 
contraction cone was designed by the method described in Ref. 7 
and was also fabricated from phenolic resin. Figure 18 is a photo- 
graph of the contraction cone and model one. The inner body was 
supported by three struts mounted in the cone as can be seen from 


the overall drawing of the assembly in Figure 19. 
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Static pressure taps were placed at intervals of two inches 
with each station having three taps one hundred and twenty degrees 
apart on the inner and outer bodies. Two stations one hundred and 
twenty degrees apart were also provided one half inch from the in- 
let plane to made radial surveys. A survey station was also lo- 
cated 8.25 inches upstream from the diffuser inlet station to make 
cross scctional velocity provile measurements and detailed measure- 
ments of the boundary layer. Hot wire ancmomecter techniques were 
used for thesc measurements and the static pressures were read 
from a manometer board. A United Sensor cobra probe was used to 
measure the inlet dynamic head and it is shown with the hot wire 
and model in " overall vicw of the assembly in Figure als 

The flow delivery system consisted of an axial compressor, two 
settling chambers, flow straighteners, and an eight inch pipe. 

The compressor was an Allis-Chalmers twelve stage axial compres- 
sor which operated at a pressure ratio of three to one. After 
leaving the compressor, the air was cooled and sent through flow 
straighteners into a large plenum. After passing through a flow 
straightener the air passed through a sharp edgcd-orifice for flow 
measurements. . Before being Supplied to the eight inch pipe, the 
air passed through the second plenum. After leaving the eight 
inch pipe the flow entered the contraction cone and then the model. 
Meriam micromanometers indicated the pressure readings for flow 


rate calculations. 
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B. EXPERIMENTAL RESULTS 

Any irregularities in the inlet velocity progile to a diffuser 
can cause separation and poor performance, therefore it was neces- 
sary to insure that the inlet ues was uniform. The hot wire 
was used to make a survey at the measuring station on the contrac- 
tion cone to measure the profile and boundary layer thickness. 
The contraction device provided a uniform inlet flow, as can be 
seen by the measurements of eRe else distributions at the sta- 
tion 8.25 inches upstream of the diffuser inlet presented in 
Figure 21. The boundary layer thickness was about .2 inches and 
the displacement thickness was calculated to be .051 inches and 
the momentum thickness .014 inches. The flow at the inlet plane 
of the diffuser is also important. It was checked with a hot wire 
and the velocit rofiles for both models were found to he rather 
mate as Shown in Figure 22. The assurance of uniform inlet flow 
was important for the possibility of Mn tas this data with the 
other existing data. Uncontrolled non-uniformities resulted in 
flow separation at the diffuser inlet in previous tests with a 
badly designed contraction device. 

The most important result of a diffuser test is the pressure 
recovery of the pedente The static taps allowed a measurement of 
pressures around the inner and outer bodies and the cobra probe 


was used to measure the inlet dynamic head. The pressure coef- 


ficient is defined as 


Cp “> CP,- PD (¢, 


where . is taken as atmospheric pressure. Both diffusers had. a 
ideal recovery doefficient of .989. ‘the non-dimensional staiit 
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pressure distribution for model one is shown in Figure 23. The 
recovery and losses are also shown on the graph. Model two also 
had good recovery as shown in Figure 24. The tests are presented 
in Figures L4 and 15. It can be seen that they agree with the 


published data presented in these figures. 
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Figure 17 Inner body and model two 





pigure 18 Contraction cone and model one 


61 





ne Sieeese ions 








a =<” 





qusudtnbsa go ydeizo.0ug 


Og einzty 





63 





lz aunols 
XDWA/A 
6 2 st € 





S 
(Ul) [JDM Wo4y 4SIP=G | 


i nd. 9 


JW4ONd ALIDOITSA 3NOD NOILOVYLNOD 





64 





s = 


Ons NI 





6G BUND eat 
7 A: : 





¢ tiaaqow 
me lis Oma AeklIO 135A 






O'} 





be JUNOIS 


01} 6° 9° b > 


1a 


67 





ec 633 

SO 9WOS SsSjoqgwhs 
9)2° =49 JONVWANYOASYSd OML ISBAQONWN 

aot X cl} =uy 


seas/wq| 2g }=W 





VIII. CONCLUSIONS AND RECOMMENDATIONS 


Diffuser performance was investigated in order to obtain a 
Simple correlation to use for design purposes which related dif- 
fuser geometry to performance. From analysis of the available 
experimental data a geometric parameter suggested by Vavra, 

_fL, (Section II, p. 17), was used to relate the important per- 
Formance quantities to geometry in a simple manner. A correlation 
relating the available experimental conical and annular diffuser 
data was finally established between diffuser geometry and per- 
formance without, considering the status of the inlet boundary 
layer. The method of utilization of the correlation in minimizing 
the diffuser length for a specified pressure recovery, and maxi- 
mizing recovery for an availiable length, was also indicated. 

Iwo diffusers were designed and tested to provide additional 
experimental data and results were found to agree with the 
developed correlation. More diffusers should be designed and 
fested in the region of the correlation where the data are not 
plentiful. The influence of the inlet boundary layer and Rcynolds 
number must be investigated, since it is felt that these two fac- 


tors are among the causes of the scatter in the derived correlation. 
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APPENDIX A 


Program Design was formulated to aid in the design of the 
diffusers mentioned in the experimental section. The purpose of 
the program was to vary the outer wall angle from sixteen to 
twenty-two degrees in increments of one quarter while always main- 
taining a specified area ratio. This procedure was done for 
three different area ratios. The program calculated omega for 
each wall angle combination by the use of a numerical integration 


scheme which employed an external function. 


List of Symbols 


me 6 C«Cinmiet area 

A2 - autlet area 

A. - area 

c - circumference 

I - logical control variable 


RH1L - inlet inner radius 

RT1 - inlet outer radius 

AHR - inner wall angle in radians 

ATR - outer wall angle in radians 

RT2 - pile’ outer radius 

RH2 - outlet inner radius 

XU - upper limit on omega integration 
XL - lower limit on omega integration 
AS - (Al1/A) *¥*3 

A6 - 2*(C/AI) 


Y - omega 


Do 





Parameter Study Programs 


There were three programs to complete the parameter study, and 
the program of each diffuser type is basically the same. The same 
integration scheme was used to calculate omega as the pertinent 
parameters were varied. PARSTUDY 1 was for annular diffusers, 
PARSTUDY2 for two dimensional and PARSTUDY3 for conical diffusers. 
The main device used was a set of nested do loops and the variables 
in these programs were identical to those in the previous programs 
except as follows. 

PARS TUDY1 
ALFAH - inner wall angle 


ALFAT - outer wall angle 
LDELR - dimensionless length 


TAT - tangent of outer wall angle 

TAH - tangent of inner wall angle 

DELP - difference between inner and outer inlet radius 
Peake LUDY 2 

LWlL - non dimensional length 

Bl ~ inlet dimension 

W1 - inlet dimension 

PARS TUDY3 

R - inlet radius 


LR - non dimensional length 


Program DATA was the most extensively used program in the 
thesis and allowed the determination of all necessary numerical 
parameters for the annular diffuser data of Sovran and Klomp. 
The numerical integration scheme was used in this program to 


calculate omega, and the same symbols in the external function 
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were used as in DESIGN. 


parameters listed in the table of symbols were calculated. 


List of Symbols 


FAM - 
AT - 
AH - 
RRAT - 


DIML - 


CPI - 


diffuser wall angle combination (family of Sovran and Klomp) 


outer wall angle 


inner wall angle 


inlet radius ratio 


dimensionless length 


area ratio 


The data was read from data cards and the 


ideal recovery coefficient 


pressure recovery coefficient 


logical control variable 


fl ee 
inlet outer radius 
inlet inner radius 
area gradient 
aiffuser length 
omega 


mean radius 


inlet area. 


ie 





DESIGN 


af 
— ez y 
-_ ‘ed _ = = 
~ Ta #2 i eo 
es << = ES N - 
= = — = . 
<= ~~ — - ats a, 
eee tY%. - ~ a -_ 
— Co co ~ o—~b 
XY YF Ke Fe «x cre 
— -— <T <= NS ie 
= Bs ze a “i oe 
Gs |. lon 
oe t j uN 
[i = = = Le Neat 
Pe = = tee 
-eSs 4 S Tyee 
i : Le — ~ 
i —™ + +> ae OG —— REL Te 
wy + oO - =o ae - . 
a os. | eee rat cy i ~ 
tane — “ee = om Ul. es 
be _— — <= << - + bes 
=D << ~~ em er ~~ me 
= — ~~ > a >< It <= 
er ~— bm SF OOH Ww ae 
wy Mo << pe RR oe mL -—_ 
® we - - <= <c ONE Ren ce 
+ C < ~ pe fe ~ @eD wu Ele 
Coe < _ b— ” + eS {>< ae, — — 
— = st AN i = + r— CeX CC ~ 
<a ec * ae eo) ee ee he eer OX Fr mele i Ce e+ 
—_—— Oe ae ee Se SI Le — Gren 
NY wer NS ee mm CS bf FP FAO DD em Ce < ot 4% 
= — “Sb oh eet toa = t es 6D - eRe b— OI de 
(es ae NR Y ROS Lc) “=~ tb uPe elit ~ ~ t—— 
~~ ec COD a + oY Nee a: *# OT —O FS ie ce C ~< te 
— a < t+ «-t-— }- -— me OT ae OH RI a5 ~ .* pe ia es — ar FT 
- Ce ie. mm om ON ee od _ Oe me Te ta ow ow ° “W) b— Sera, 
Cc ee WRT mO™N YS em LK a a Ose ion C3 oe el MN —~ Fa =) co kl 
~— ~ eo IL EeN em | eUrmmO mL, a = 3 Pec —~ 6 30m che to ¥ 
RW DW ewe~w aS oN oe + See Gn oUt =e tfh,.*e mnie Ol © MUN AK OX 
s/t NS Ns Weare ZT at+ TNOFG TD WOON OM OC HC et em AN FN 8 6m a ee + 
< HOOD ee Neem FOO YN NL Om OOO LNA A Owe CN Oe COO NN CO mt SO et 
wa: ae el OF - MWIO I HAH Lee —= — | TMM BOCAS Om + owl Lil © © © — @ e-- _. 
eH me OS i be de He ie i ee ee er ew OR Cel es QQ OOS be Oe OD. 
hd od HUSH eet ree HOW eH weeconxeD + ft ot RR TNS TR Fe te OIE WHhrkha~~— 
RST Tt I RR HONE NL towel) » @ l>O> Neo Sl —— See H ' aac Zea tt 
>< us Wee TNS ee a < il HOCH HCO Dee Her DOH Oi HO ft ON He VQORrFLD NSE CULuUnONCrZ Se Dene 
Wi Oar a IKE CCU OI OAC KE OSM KD Ue OO OS a LO Ee ee DON UL Ot ox NK 
<< can en 
aad < 
WY iL 
mei OY c aoa 


n; Sue a DM OG 


1 
C 
< 


72 





SSS NG) 
~N SN 


sf 


F760 2S 


CONTOAL “Pal ire R STUDY 
PARS THUD 


<4 


L¥4(L) 


ae 


—e (|e € Sainces 
mo 
eo~JN YD 


ee > ep SS 
mS AG) 
2S = 
o =~ 
bet peed pet 
~~ — OM 
=> 


~_, 
“we 


ZN PNA Dat 
“S@ em see tl Co df 


ao tw tt Sry 


FAR) / 7COS CALE AR) 


A e(PR+EXU TAN) €¥*2 
La/AR ED 


Rape 


jc ree ‘Tp 

Wee Tos (rCTICATC, TANT ECT CNR=C TAN) 
AYRE 15620 

Me a 2 (ECT CAFC pT AN) *+FCTCA-Cy,TAN)) 
ee ee nk a a 
ee POLE SS Ui ag ae ae 


1735" (FOTCAtC, TAM) + FCT CA _ TAN) ) 
Meet tate VOVM(As TAN) ) 

at oer as) See 

MEGh=" « Cie 02%, "LA RE Mer Oe ech Ty WEE A=" ab ee 


va 


a OP > ti 


at 
ee 
| 


<0 @ © >ce@ ff melt fl 


< @ 
Be eR es AL 


C~ >: 


o. N= oe? 


Q4Cm PAST x 


— Oi mah te em ans: 


2eHo4os5nrsre~ Sew ht WOW roe oe ee ee 


O 


—* 
ze 
ee 


<€m7aewn>< Nx SAS TT 
Naz TOWN Ihre Hh bee 
= JAH Pw tha 
TCuMe xe pe 
A ee 
VD are >We RN 
8 ee 


— 
t 


—? 
Kom 


73 





; THO DIMEMSTONAL PARAMETER STUDY 
C PAR STUDY? 


Pa 
7 


lew The oO 
=) 


a 
ae 
oc =>- 
cnr — 


Ne 
” 


ree eee OK 
eo ee fT] um 


oo 
ad 


NN eee [Pp ree TZ OS 
7} 
wm 
Es 
aoe 
“™ OJ 
2 
/ 
= 
een 
ot 
>> 
WD 
— 


bat pot tat @ Of) 

e + 

i pe Be ere re iT It it 
fre et ee pt ted pee 


r= @ 
e ar 


’ 
+e 


Bie XE) 

a0 

oo Manohla hid +PCT CASE, Ia 1 a) 
DO*CECTC ALC TAN La ee 1 kOe Wah, Bly Wi 
Beater et TAN, B1lyALT)FFCTIA-C, TAN, BlyWwl)) 


tC | iA» Ae, ae 


a. Ce cA ALR a Poe cee 


Me YP Wee mean 


<@ o e 
ee Oe 
ae 

w 


“Ms Mw 


ne me re Noe 
wom oo ee oe 


fmWenX et hPU HEMT UW ose tl PpwNrer lle >in 


Wot a ao ea a See 


gree AS USI UT yo 


on 7 


qua ONC oe 


ee 
pont a ye 
a aon 
ee 
[I 
Rev 
e 
aa 
x) 
>< 
* 
™~ 
oc 
| 


NNO Beal) See = 
~~ 
> 
e 
Ne 
IND 


—_ 
Nal 

QI HNN ZW Drea sas <e 
Cs _—l + 


co 
=> 
ili at. ade 

>> 

oe 

NQ 

2 

~ 


— ——- 
en 


m™ 
~- 


“laa < 
{v 
@ 


sonnel oye ee 

beet lt © ——t jf © mat ee | PTA 
—- 
« 


a Sheet 


os amy * 
— 


a ee 
2S (Sea Se Se 
th» 


ZH Te am Re MD 
SN TS Fe —t 


mH ITKOM DPT 
Preeti th eee: [re 
Jam NUN PW. 


—~, 
reo 4 


//GOD. SYSIN DN » 


74 





2> 


— 


—) 


_— 


femwed 


ANNULAR PARAMETER STUDY 


2))/A) 


CT CATE gS 1] sie ls TA OC aC ls ee 


KS) + X25%, 


PARSTUDY3 
Miri Giy REAL*4(1 } 
EXTEVPNAL FCT 
WRITE(5,25) 
FORMAT( #16) 
ere oll. 
Ree Ald =i fs 
(ee= on / 57, 3 
eae =).. 
R= Ave 
ioe iw | 5 
NG 2 eat 14 5 
N=J+1 
Pits SEN. 2 Py 
DN 4 [=)].4 
pea OL / Beet [eZ ye, 5 
NELREPO—P] 
AL FATR=ALFAT*5 7.3 
ALE AHR=EALEAHS5S 7.2 
Pat tC aLe AT) /COSCALEAT ) 
TAHHESTIHCALT AH) /COS CALF AH) 
XUSH=LNLLR*DEL®P 
AP=H(OBVGL*( (2M%+ XU TAT) 2 2-(R [+ XUX TAM) KX 
Oe — (laff? ) ¥¥ 2 
Pre =2. eRe TAL ole? Ps TAY 
Piet ( (TAT )#*2 = CT AH aes Peo, 
les pend 
A= eS CAUTAL } 
H=KI- Xl. 
C= 2484 199) *° 
Meee TLR CECT IASC, SRI, RU, TACT (A-C,RI,RO, 147. 
hy 
Cee) 150*1 
ey 4G peer Oe, TAT, TAR) FFOTCA-E€,R81,°99,TA 
Wy Tot) ) 
C=.72 66R57%D 
VYsYVte L326 53% (FE CT( AC, RI],RU, TAT, TAH) FCF CTI(A-C,RI RO, TAT 
Peete} 
GewiGici cod a 
ee b>] PSs (Fr 
eA) 
Wepe{ Vt. 1651197 CICA 21y,RO,TAT, TAY) } 
tna ibts . YAN 
SNRAAT CC? ©, \O*MEGAS', FSe2) 3X, *LDELR=!,F 6029 3X; *DADX= 
So.) ; 
Pale m1 +7. i) 
CONTINUE 
Pal = ea ; 
mig | FAT+5./ 5764 3 
eo NUE 
Lae FUNCTPIN FC 
ALFAHSALEAH+5., /57. 3 PUNT LLIN 
ALFATHALFAHF+1Le/ 5702 AL=?2. 6 
CONTINUE ete en 
ee c oe ee 
ALEAT=16°/5743 Ke ee 
LOELR=LDELR +2. Kaa ce* 9 14] 
CONTINUE (THY xs 
€7 7p . Ne Ua 
pares ee 
Sa S¥STh py 


9) ) 





ANO KLO‘“P 


SOVRAN 
MY) FE CTCA-C gah Re Aneta 


+ 
\ 


FROP 
wARTLISX=(F CTC AFC, ATR,AHR EF 


MEDUC ie 
FAT Ne 


lS Cee) 
NS CAH® 


DATA 
? Ot: XK ’ : NEW 
GO Tus S 


GN) TO 4 
FCT (X cATR AFR ee 


a eS 
OVS sl CCL Na 
oe TOGA <T we CO en 2 


FO DTA Cw ew Vw OSL 


) 


Ne NOOO «0 ect<in it 
teccuod | RNS ~~ 
AMPH erUWNLI Ls 3 


. * 
LO AO ee mee pen pl bee SS st ert sun 4. UAW w-t 


Com 
— T- 
am 
==" © €¢ © 
VNC OG an ) 
eC:C 3t te 
Nits se > 
oe ee 
a2 GF aM 
mie << 
Co | 
fmm EU N 
ee es OK 
Somes + 4 
may 


4 


—' @es ff Ij) b- o 
be be oO a OX 
ofl cn 


Came fm fl me ee Oe FEL! 
Set et Soe oe TT 

. -.- * -_ 
er DmAAINAY NCR. Tee an et rae 


ete ye te CO a t. 
“LEMS SOL VSTi Nh oh 
eee OS Co OS Se 


tet FQOU3B2SCR (FCT (At+C ATR, ARR RHI IFFCI(A-C AIR, Smee 


4) 40155) 
0 3 668574 
te) 3 3°538 (CF CTEACL ATR. AHR. RHI) FFOTOA-C.ATR «AHR. PH 


Bee yueXL) 
4 Ra’ ALLEN 


ere eee tt 
DR REAM 
Lo x x™mM< 


fe ett TOY ILI YR TO 
Se lie Ge i te) ee er SO et ee, df 
We OY STUY OCA ORK PEE OOK KS 


D=XU-XL 


ty seCyemu NOOK ONY tea we OY weer 
, > OS CCS SUL OMY OU DO DUES 
Aan gRaawes eowbK CEU RH OnwuU 


99 
Tai. 


76 





sei mis . . ° be ee 
Co Tr — “ix aso a be et. mo. she em 


ont ‘. = 

a ee aa eer Gk GS cn ec h 

ALOT GENRE ON OS meg ECE MPN ET OS 
ts Cit. 

e@ 


a ee ey fa et, oe ete COA INC 
SEND Me NIU re NLS OO CL OTN LT 
Mi 7 DINO eee Gamo ea 
ACR SR COMMO Sy tS SS IAI OR ENNIO SS we 


d 


SLO ONE COC Ot FOI OM MONG 


SO ere eg et tc emt ee ee te 
NOES) SOC NDNA OCA TOR MANL OAS TE OL OSCAW OD OT SERV DO EON NM OMS MOG AOR MES _ OW amen 
SR MAT MR REO ST ONMUAH DM OO Ce tS TORO URAN SAM OO mt ROP RS OE Re NOAM ST ITN™ eee 
AMA OSE OO OMEN CORR ROO MHORR ROPER SEO MRM ONS DN OR RO MOVIE OSU CRE OORK Goes OC Our oO 
roe ee ene once es ac? omen) * 8s 8 fe 6 8 we eo be eee ee eee ee oreo wee eo 
ee CA OO a Or AMANO S Corte VS NCR CNC eMC Mmm. AS 
LaOewsaraHeanahy COCm™ oe (Gwe eee he Or es OO sige < AAO ASL OOD OR aa ge: Mais: 
pe ee eo TT meee OC case RON OC OIE et OSA LAD ted mI CIF OUT OO Oe iho igs POCA CA RaN ORM SED 
oeeweenanaee5uscegeeee8teedopee#e0ecetkv#2#eegeegesv#seeoqeee#eeceererveeeee#e#ee%e*#ete#eteegegqee#t*e¢eeeteeeoeeetee#ee 8 © @ @€ O@ ¢ 


in a ee a SS ae ON NO Net ee a NO oe eee 

GO CUNS Pe Gaertn UN AS SI OG or ON em ee oT Ne coe CU nt? Sent ON ee Nike ee 

NOC m™m™ _ -Si.gc™ PAHS ORT AOMUS ER MOOIMR ES GMMR Tee Pelco ~~ CSO OAl Ole Orr Or ls TSS OL row n 

Cem ee eM Oe em ee ea me sear eo ee HereeHEHoearesertesHeroereeesoeeroeeeHics 

CON Nepal aaa: WVU ONIN TN OR SEE UNC ON ONC Moa ee tees eee as ae oe Se ae a OS es: 
Sane gman pond apt on 


’ 
a oe 2 See > *e : ~ 


ei. 5) +X Bae 6) 7K) 


Rate ree ane et ne Pe pte € (2 25 ° 2 : 5 
cS ao aes Noe \ 0 eee Ca, Se eee 5 ee eee. ae c. ie inc ee tee Ace in ic ioways mnie tA 
ee ee eee a nines Pieters WO i ati 
OS a ee cae ee le er Bee See TX. ee ee Of? et oes See es tee 
>< 0 tt C2 eee oe eC a A or ro a hee ores. = per (ae ae” Fs. = a wel a. “" > 2 Cian. ; 
—< oe se ere Se ere oe ea NNO S- a Vatattsl tates nian Ser Fee ee OM surat Co ONL VU ee» 
— > Eo er | ee Ae e > re rN rt ee FB ec wets ce eC te oe 8 
StI AINA MARAE IAL RD = INOOGC GROG MM DAI OTT SOS TEL VOT YN eS oe Ghee 2 Se TON It: 
ec eriente nto e. ra ae e es age Ge one cial eee el ee < a ce med Am ee a ait. Epes Stet ie pao ere ae .- E. rete roe- 
be eee ee eS Lr BES ee ee a oe tT OU SU 
ae a 1 ot ae oh de OE a ee Seca SC Me eee Se Oe ee In ee rr te 
ee eee ee es eee Semen) - Gea = oo gnc + tts 2 Fa + vito ue et 8 oe 8 oe 2 
ee a iG AGN ae I eS Oe ee Rs LOL LOUD _ we Bs 
D9 mt mt es ab et mt RMR AAG OU AION: NINA es: CELA OAT OSS Rea a“ me Se Cae aia oe age: RUC aes Ci. ew re 
Oli WC Ware IN ene LU tet i eee” Gi OS CO UNI LL O- aS ee ee ee se eo ae ee ee A Ct — 
CASRN sede ese ed es aS rs mm ed i et edt ft eee Ce cM ge aes ees SINS ee Ne we ON ee r ee Noss aera De - a sta an pny ey eed eed esi ot ae ae port ee et a  - 


ee ee ee ines ote risie ce tee ee hh eect c isd vo spetesn> 
Saba el ee er ee ae eM Ne ee ae ce ye EA ee LOIN See 2 eee 
Rt et et pd ed met MOU OA QS OI AIO CU AI OU AI AI AI QUCU AIOIINUAIAIONIOGN AURA OOOO DI I RR tet HH RRMA NN NN NOIN Go rm a 


v7 





BIBLIOGRAPHY 


Bocanower, A. B., Kline, Ss. J.,7and gohnston.. J. P-. 
"A Unified Method for Predicting the Performance of 
Subsonic Diffusers-of Several Geometrics," Report PD-10, 
Thermosciences Division, Mechanical Engineering Department, 
Stanford University, May 1965. 


Cockrell, D. J. and Markland, E., "A Review of Incompressible 


Patiuser Flow,” Aircraft Engineering, v. 35, pp. 288-292, 
Be ecoper 1965". 


Gapp, D. R., A Design Criterion for Arbitrary Diffusers, 
Aeronautical Engineer Thesis, U.S. Naval Postgraduate 


School, Monterey, California, June 1970 


Gleason, J. G., An Investigation of Current Subsonic Diffuser 
Desi Knowledge with Emphasis on Jct Engine Annular 
Diffusers, A Pratt and Whitney Summer Faculty Program 
Project, Summer 1963. 


Remeau, L. R., Johnston, J. P., and Kline, S. J., 
"Performance and Design of Straight, Two Dimensional 
Diffusers," Report PD-8&, Thermosciences Division. 
Mechanical Engineering Depavriment, Stanford university, 
September 1964. 


Sovran, G. and Klomp, D. D., Experimentally Determined 


Optimum Geometries for Rectilincar Diffusers with 
Rectangular Conical or Annular Cross-Section, 


Engincering Development Department Research Laboratories, 
Gencral Motors Corporation, Warren, Michigan, 
16 November 1965. 3 





Tsien, H. "On the Design of a Contraction Cone for a 
Wind Tunnel, "Journal of Aeronautical Scienccs, 
pp. 68-70, February 1943. 


Wawra, M. H., Agard Lecture Series 39,"Lleepure Series on 


Basic Elements for Advaneed Designs of Radial-Flow 
Compressors, June 1970. 


78 





un 


INITIAL DISTRIBUTION LIST 


No. 


Defense Doeumentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Chairman, Department of Aeronauties 
Naval Postgraduate Sehool 
Monterey, California 93940 


Kyriaeos Papailiou (thesis advisor) 
Dip. Mech. Eng., Dip. Aero. 

von Karman Institute 

72, Chaussee De Waterloo 
Rhode-Saint-Genese 

Belgium 


refessor M. H. Vavra 
Department of Aeronautics 
Naval Fostgraduate Sehool 
Monterey, California 93940 


K. W. Sharer 


6015 Lamont Court 
Springfield, Virginia 22152 


7S 


Copies 








Security Classification 


DOCUMENT CONTROL DATA- R&D 


(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classttiec) 


ORIGINATING ACTivITY (Corporate author) 


28. REF ORT SECURITY CLASSIFICATIO?*s 


| Unclassified 
Naval Postgraduate School 26. GROUP 
Monterey, California 93940 
3 REPORT TITLE 
An Investigation of Subsonic Diffuser Performance 
ja OESCRIPTIVE NOTES (Type of report arid, inclusive dates) “ 
| 
S. AUTHOR(S) (First name, middle initial, last name) 
Kevin Woods Sharer 
. REPORT OATE : - 7a. TOTAL NO. OF PAGES : 7b. NO. OF REFS 
June 1971 81 8 
Pea. CONTRACT OR GRANT NO. Toa. ORIGINATOR'S REPORT NUMBER(S) 
b&b. PROJECT NO. 
C:. 9b. OTHER REPORT NOS) (Any other numbers thet may be assigned 
this report) 
d. 


eee 


10. OCISTRIBUTION STATEMENT 


Approved for public release; distribution unlimited. 


| 11. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


£13. ABSTRACT 


This report investigates subsonic diffuser performance, with emphasis 
on conical and annular geometries. A correlation is presented which ajds in 
the prediction of performance. ‘Two annular diffusers were designed and 
tested to help substantiate the correlation. 


DD FO 4473 ~~ (Pace 1) 


S/N 0101-807-6811 80 ri “Secunty Classifica) 





Security Classification 


7 
| ae 


eta 


LU 
ular diffuser 

fuser 

fuser performance 

e parameter 

ical diffuser 

elation for diffusers 






Li (FORM 4 473 (BACK } — . 


H101-507-682! 81 Security Classification Wes 2 




















Thesis 1276406 
$4345 Sharer 
Con An investigation of 


subsonic diffuser per- 
formance. 





